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ABSTRACT 

We present a detailed study of the long-timescale X-ray variability of the Narrow Line 
Seyfert 1 Galaxy (NLS1) MCG-6-30-15, based on eight years of frequent monitoring 
observations with the Rossi X-ray Timing Explorer. When combined with the pub¬ 
lished short timescale XMM-Newton observations, we derive the powerspectral density 
(PSD) covering 6 decades of frequency from ~ 10 -8 to ~ 10 _2 Hz. As with NGC4051, 
another NLS1, we find that the PSD of MCG-6-30-15 is a close analogue of the PSD of 
a Galactic Black Hole X-ray binary system (GBH) in a ‘high’ rather than a ‘low’ state. 
As with NGC4051 and the GBH Cygnus X-l in its high state, a smoothly bending 
model is a better fit to the PSD of MCG-6-30-15, giving a derived break frequency of 
7.6I3 0 x 10~ 5 Hz. Assuming linear scaling of break frequency with black hole mass, 
we estimate the black hole mass in MCG-6-30-15 to be ~ 2.91R x 1O 6 M 0 . 

Although, in the X-ray band, it is one of the best observed Seyfert galaxies, there 
has as yet been no accurate determination of the mass of the black hole in MCG-6- 
30-15. Here we present a mass determination using the velocity dispersion (Mbh — cr*) 
technique and compare it with estimates based on the width of the Ha line. Depending 
on the calibration relationship assumed for the Mbh — o’* relationship, we derive a 
mass between 3.6 and 6 xlO 6 Af 0 , consistent with the mass derived from the PSD. 

Using the newly derived m ass and break t imesc ale, and revised reverbera¬ 
tion masses for other AGN from iPeterson et alJ (120Q4T) . we update the black hole 
mass/break timescale diagram. The observations are still generally consistent with 
narrow line Seyfert 1 galaxies having shorter break timescales, for a given mass, than 
broad line AGN, probably reflecting a higher accretion rate. However the revised, 
generally higher, masses (but unchanged break timescales) are also consistent with 
perhaps all of the X-ray bright AGN studied so far being high state objects. This re¬ 
sult may simply be a selection effect, based on their selection from high-flux X-ray all 
sky catalogues, and their consequent typically high X-ray/radio ratios, which indicate 
high state systems. 

Key words: black hole physics - galaxies:active - galaxies:individual:MCG-6-30-15 
- X-rays:binaries - X-rays:galaxies 


1 INTRODUCTION 

It is now reasonably well established that the X-ray pow¬ 
erspectral densities (PSDs) of AGN are broadly similar to 


those of galactic 

black hole X-ray binary systems (GBHs) 

iM c Hardtll98& 

Edelson & Nandra 1999: jlJttlev et al. 2002; 

iMarkowitz et al 

12003; M c Hardv et al. 2004). The PSDs are 


described by powerlaws of the form P(v) oc v a , where P(v) 
is the power at frequency u, but a varies with frequency 
(see Section 3.2.1 for details). At high frequencies the PSDs 
of both AGN and GBHs are steep (slope, a, ~ —2) but be¬ 
low a break frequency, vb, typically at a few Hz for GBHs, 
they flatten to a slope of a ~ —1. To first order, the break 


timescale scales linearly with the black hole mass but there 
is still considerable uncertainty in the mass/break timescale 
relationship and hence in our understanding of the physical 
similarities between AGN and GBHs. 

The uncertainty arises largel y because G BHs occur 
in se veral distinct states (e.g. see iMcGlintock fc Remillardj 
l2003h . Most commonly they are found in the so-called ‘low’ 
state where their X-ray fluxes are low and their medium en¬ 
ergy (2-10 keV) X-ray spectra are hard. The second most 
common state is the ‘high’ state where their fluxes are high 
and their 2-10 keV spectra are soft. The PSDs of GBHs in 
these two states are quite different. I 11 the low state, there 
is a second break, about a decade below the high frequency 
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break llNowak et al.ll99sib Below the lower frequency break 
the PSD flattens further to a slope of zero. In the high state 
there is no second break and th e PS D continues with slope 
~ — 1 to very low frequencies llCui et al ] 1997)- In addition 
the break from slope ~ —2 to slope ~ — 1 occurs at a higher 
frequency in the high state than in the low state (~ 15Hz cf 
~ 3Hz in the best studied GBH Cyg X-l). Before comparing 
AGN with GBHs it is therefore important to know whether 
the AGN is a high or low state system. 

Based on their 2-10 keV X-ray spectra, it has gener¬ 
ally been assumed that AGN are the equivalent of low state 
GBHs. However, although the high frequency {v > 10 _6 Hz) 
parts of many AGN PSDs are now reasonably well deter¬ 
mined, the lower frequencies are, in general, not well deter¬ 
mined and so it is not usually possible to be sure whether 
AGN are the analogues of low or high state GBHs. In gen¬ 
eral it is not possible to determine whether there is a second, 
lower frequency, break or not. The difficulty in determining 
the low frequency shape of AGN PSDs has been in obtain¬ 
ing well sampled lightcurves stretching over sufficiently long 
timescales. For an assumed linear scaling of break timescale 
with mass, we require well sampled lightcurves of a few years 
duration to detect the second, lower, break in a black hole 
of mass ~ 10 6 Mg. 

Prior to the launch of RXTE in November 1995 it 
was not possible to obtain long timescal e lightcurves of 
suffic i ent quality but, since early 1996 we llM c Hardv et alJ 
[l998|; |uttley_et_alj |2002|; lM|!Hardy_eh-a]J an d oth¬ 

ers (IEdelsoi^r^andralll999nMarkowit^e^ni2003ll . have 
been monitoring a small sample of AGN, typically observ¬ 
ing each AGN once every 2 days, and so are able to deter¬ 
mine the shape of the low fre quency PSD to high p recision. 
Using these observations we dM c Hardv et aljEfiolIf recently 
showed that the PSD of the narrow line Seyfert 1 (NLS1) 
galaxy NGC 4051 was, in fact, identical to that of a high, 
rather than low state GBH and so provided the first definite 
confirmation of a high state AGN. [Markowitz et alJ • 20ll.'»i 
plotted the black hole mass against the timescales associated 
with well measured PSD breaks, some from slope ~ — 2 to 
~ — 1 and some from slope ~ — 1 to ~ 0. They claimed that 
there was a linear relationship between the two parameters . 
With a somewhat larger sample we dM c Hardy et alJ 1200 4') 
plotted the black hole mass against the timescale associated 
with the specific PSD break from slope ~ — 2 to ~ — 1 for a 
sample of AGN. We found that the best fit to all of the AGN 
together had a slope flatter than unity and did not extrap¬ 
olate to either the high or low state break timescales of Cyg 
X-l. However we noted that the break timescales associated 
with broad line Seyfert 1 galaxies were, for a given black hole 
mass, generally longer than those for those objects usually 
classed as NLSls. NLSls are probably not a distinct class of 
object but, more likely, just lie at one end of a spectrum of 
AGN properties, characterised perhaps by a higher accretion 
rate. However taking the crude broad/narrow line distinc¬ 
tion which is commonly used in the literature, we note that a 
linear BH-mass/break-timescale relationship fits broad line 
AGN and Cyg X-l in the low state, and a displaced lin¬ 
ear relationship fits the NLSls and Cyg X-l in the high 
state. We therefore suggested that there is no fixed break 
timescale/mass relationship which fits all AGN but that the 
break timescale/mass relationship may actually vary in re¬ 


sponse to one or more other factors, which might be accre¬ 
tion rate or black hole spin. 

In order to test and clarify the above hypothesis, it is 
necessary to place AGN on the break timescale/mass plane 
with high precision and to determine whether they are high 
or low state analogues. There are very few AGN for which 
the X-ray observations are extensive enough to make the lat¬ 
ter determination possible and a particularly important ob¬ 
ject is therefore the Seyfert galaxy MCG-6-30-15. Although 
MCG-6-30-15 is not always described as an NLS1, it has 
many of the properties commonly associated with NLSls, 
ie rapid X-ray variability and relatively narrow permi tted 
emission lines (FWHM 1700 km s -1 dPineda et alllOSllf c.f. 
the 2000 km s~ 4 limit which is usually quoted for NLSls). 
In this paper we therefore class MCG-6-30-15 as a NLS1, 
like NGC4051. 

MCG-6-30-15 is one of the best studied X-ray bright 
Seyfert galaxies. It was the galaxy in which the first detec¬ 
tion of a relativistically broadened X-ray iron line was made, 
providing di rect ev idence for the existence of a massive black 
hole llTanaka et al.lfl995t) . We have monitored it extensively 
with the Rossi X-ray Timing Exp lorer ( RX TE) an d have 
presented its long timescale PSD (lllttlev et ai]l2003b based 
on 3 years of o bservations (1996-1999). From that early, lim¬ 
ited, dataset llJttlev et alJ ( 2002 4 were able to show that a 
simple unbroken powerlaw was not a good fit to the PSD 
(at greater than 99 percent rejection confidence) but that 
a broken powerlaw did provide a good fit (rejection con¬ 
fidence 33 percent). However it was not possible to deter¬ 
mine, with any accuracy, the low frequency PSD slope. If 
they assumed a powerlaw of slope -1, Uttley et al. were able 
to derive a break frequency in the PSD at ~few xlCU 5 Hz 
and to estimate the high frequency powerlaw slope (~ —2). 
Thus Uttley et al. were unable to distinguish a high from a 
low state system, or measure the break timescale with high 
precision, although they preferred the high state interpreta¬ 
tion as a low state inter pretation implie d a lum inosity close 
to the Eddington limit. iVaughan et al. ( [2003 1 have subse¬ 
quently published a detailed study of the short timescale 
X-ray variability of MCG-6-30-15 using observations from 
XMM-Newton. These observations are insufficient, on their 
own, to determine the low frequency PSD slope but, by also 
assuming a low frequency PSD slope of -1, they have deter¬ 
mined the break frequency to be ~ 1 x 10~ 4 Hz. 

In this paper (Section 0 we present our full RXTE 
observations, of 8 years duration, and consi sting o f over 
800 se parate observations, compared to 100 in llJttlev et alJ 
d2002dl . In Section ED we combine these RXTE observations 
with the published XMM-Newton observations and derive a 
PSD of excellent quality covering over 6 decades of frequency 
from < 10~ 8 to > 10~ 2 Hz. We determine accurately the 
shape of the PSD below the ~ 10 _4 Hz break and we thereby 
show that MCG-6-30-15 is, like NGC4051 (lM c Hardy et aid 
l2004lh the analogue of a high state GBH. Using our improved 
determination of the low frequency PSD shape, we slightly 
refine the break timescale. 

The most widely accepted method of determining black 
hole mas ses i n Seyfert galaxies is that of reverberation map¬ 
ping fe.g. lPetersonl200Jl but this technique has not yet been 
applied to MCG-6-30-15. The best current estimate of the 
BH mass in MCG-6-30-15 is ~ 1 x 10 6 Mq dlJttlev et alJ 
l2002lb 
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Figure 1. Long Term RXTE 2-10 keV lightcurve of MCG-6-30-15. Each data point represents an observation of ~ 1 ks. 


That m ass is based on an estimated bulge mass of 
3 x 1O 9 M 0 llRevnoldj [20001 and on the correlation be¬ 
tween black hole mass and galactic bulge mass presented 
by EE n del ( 1994) . where he claims that Seyfert galaxies 
have a lower black hole mass, for a given bulg e mass, t han 
was cl aimed in the ori ginal relationsh ip of lMagorrian et ahl 
i 1998h . In a later paper WandeJ ll2002il revises the black hole 
mass/bulge ratio to 0.0015 and shows that Seyfert galaxies 
fit the same relationship as normal galaxies. Using th e re¬ 
vised ratio, and the estimated bulge mass from iRevnolddl 
feOOOlf . the revised black hole mass would be 4.5 x 1O 6 M 0 . 

An upper l imit to the black hole m ass of lO 7 M 0 has 
been derived bv lMorales fc Fabiani (‘ 2002;’) . using a method 
based on balancing the radiative and gravitational forces 
acting on outflowing warm absorber clouds. However it has 
recently become clear that the successful technique of deter¬ 
mining black hole masses from the central velocity dispersion 
of the galaxy bulge can be applied to active, as well as qui¬ 
escent, galaxies by observing the width of the Call triplet 
in the far red, where the AGN continuum is not dominant 
llFerrarese et al • 2001). In Section 2] we report the deter¬ 
mination of the black hole mass in MCG-6-30-15 from cen¬ 
tral velocity dispersion measurements. For comparison (Sec¬ 
tion EJ we also estimate the BH mass from the width of the 
H/3 emission line and the con ti nuum fl ux a t 5100A, using 
empirical relations from liEsni et all f 20001 ). With greater 
error, we also estimate (Section EJ the BH mass using the 
‘photoionisation’ method ( IWandel et a DIEII- 

In Section |7| we discuss the implications of our results 
for the comparison of AGN and GBHs. 


2 X-RAY OBSERVATIONS AND DATA 
REDUCTION 

2.1 RXTE 

The RXTE observations discussed here consist of our own 
monitoring observations, which have continued since 1996, 
and two ‘long looks’ in August 1997 and July 1999 which 
we have obtained from the RXTE archive. The monitoring 
observations cover timescales from less than a day to ~few 
years and the long looks cover timescales from ~days to 
^minutes. 

2.1.1 Monitoring Observations 

The RXTE monitoring observations, of typically ~1 ks dura¬ 
tion, were made and analysed in exactly the same way as for 
NGC4051 <|m! Hardy et al ]». Prior to 2000, we used a 
quasi-logarithmic sampling pattern, covering all timescales, 
but in order to improve the S/N on our resulting PSDs, we 
then increased our coverage. From 2000 onwards, our typ¬ 
ical observation frequency has been once every two days. 
In addition, to sample properly the higher frequencies, we 
observed once every 6 hours for 2 months in 2000. Each ob¬ 
servation, of ~ lksec duration, contributes one flux point to 
the monitoring lightcurve. We do not attempt to split the 
monitoring data into higher time resolution bins. 

The observations were made with the Proportional 
Counter Array (PCA) which consists of 5 Xenon-filled pro¬ 
portional counter units (PCUs). We extracted the Xenon 1 
(top) layer data from all PCUs that were switched on during 
the observation as layer 1 provides the highest S/N for pho¬ 
tons in the energy range 2-20 keV where the flux from AGN 
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Figure 2. The 2-10 keV RXTE lightcurve of MCG-6-30-15 cov¬ 
ering the two-month period of four observations per day. Each 
data point again represents an observation of ~ 1 ks. 


is strongest. We used FTOOLS v4.2 for the reduction of 
the PCA data. We used standard ‘good time’ data selection 
criteria, i.e. target elevation > 10°, pointing offset < 0.01°, 
time since passage through the South Atlantic Anomaly of 
> 30min and standard threshold for electron contamina¬ 
tion. We calculated the model background in the PCA with 
the tool PCABACKEST v2.1 using the L7 model for faint 
sources. PCA response matrices were calculated individu¬ 
ally for each observation using PCARSP V2.37, taking into 
account temporal variations of the detector gain and the 
changing number of detectors used. Fluxes in the 2-10 keV 
band were then determined using XSPEC, fitting a simple 
powerlaw with variable slope but with absorption fixed at 
the Galactic level of 4.06 x 10 20 cm -2 dElvis et, alJll98flll . 
The errors in the flux are scaled directly from the observed 
errors in the measured count rate. 

As with NGC4051, we produce a lightcurve in flux units, 
rather than raw observed counts/sec, so that we may use 
together data from periods when the number, and gain, of 
the PCUs changes. The resulting lightcurve, from 1996 to 
2004, is presented in Fig. Q 

From the monitoring observations we make three 
lightcurves. The 6-hr sampled lightcurve consists of the 2 
months of observations 4 times per day. The 2-d sampled 
lightcurve consists of the period from 2000 onwards where 
sampling is once every 2 days. The total lightcurve consists 
of all of the observations since 1996, but when used in the 
PSD analysis it is always heavily binned to 28 or 56-d reso¬ 
lution. 


2.1.2 Long Looks 

There have been two long looks with RXTE , one in August 
1997 and the other in July 1999. The total duration of each 
observation was about 9 days, consisting of continuous seg¬ 
ments of ~ 3ksec, separated by gaps of approximately equal 
length caused by earth occultation. The August 1997 obser- 



0 2xl0 5 4x 10 5 6xl0 5 

Time (secs) 


Figure 3. RXTE background-subtracted lightcurve of MCG-6- 
30-15 in the 2-10 keV energy band, with 128s time bins. The 
lightcurves are shown as histogram segments. Gaps are left where 
there are no data. Bottom panel is the August 1997 observation 
and top panel is the July 1999 observation. In both cases the time 
is measured since the start of the observation. 


vations have been published bv lLee et al. ( 2000 ) who have 
carried out a pre liminary analysis of their variability prop¬ 
erties. iLee et all d2000l) clai m a break in the powerspectrum 
at ~ 4 — 5 x 10 -6 Hz. iNowak fc Chianel d200(t) have anal¬ 
ysed the same dataset, together with ASCA observations, 
and claim that the resultant PSD resembles that of a GBH 
in a low state. The claim that the PSD is flat below 10 _5 Hz, 
has a slope of approximately -1 from 10 _s Hz to ~ 10~ 4 Hz 
and, above ~ 10 -4 Hz it steepens to a slope of -2. Neither 
Lee et al. or Nowak and Chiang ca rry out sim ul ations to 
determine the validity of their PSDs. lUttlev et al] (l2002h do 
carry out simulations and conclude that that although there 
is a break in the PSD of MCG-6-30-15, the model of Nowak 
and Chiang can be rejected at 99 percent confidence and 
that there is significant power in the PSD below 10 _5 Hz. 
A time-series analysis of the July 1999 observations has not 
yet been published. 

The two long looks were analysed in exactly the same 
way as the monitoring observations, except that we retained 
16s time resolution. As the July 1999 lightcurve has not 
previously been published, we include it here (Fig. [3J and, 
for the convenience of the reader, we also include the August 
1997 lightcurve which has previously been published by Lee 
et al. 
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Figure 4. XMM-Newton background-subtracted lightcurve of 
MCG-6-30-15 in the 0.1-10 keV energy band, with 20s time bins, 
using data from the PN CCDs. 


2.2 XMM-Newton 

The XMM-Newton timing observat ions of MCG- 6 -30-15 
have been discussed extensively by IVaughan et alJ (|2003h 
to which we refer readers for a detailed discussion. For the 
convenience of readers, we reproduce their lightcurve here 

(Fig®. 


3 PSD DETERMINATION 



Frequency (Hz) 


Figure 5. Long timescale RXTE PSD of MCG—6-30-15. The PSD 
is well fitted (P=68 percent) by a simple powerlaw of slope 0.9 it 
0.25 over three decades. The underlying undistorted model PSD 
is shown by the smooth continuous line. The distorted model for 
each of the three datasets, together with its errors, is shown by 
the individual points. The observed dirty PSD is given by the 
jagged lines. 


3.2 Combined Long and Short Timescale PSD 

3.2.1 RXTE monitoring observations and XMM-Newton 
4-10 keV observations 


The method used to determine the PSD is the same Monte 
Carlo simulation-based modelling technique llUttlev et alJ 
l200di) . PSRESP, which we employed in the analysis of 
the combi ned XMM-Newton an d RXTE observations of 
NGC4051 JM C Hardy et alJl2004l) . This method is able to 
take account of non-uniform sampling and gaps in the data. 

3.1 Long Timescale PSD 

A variety of datasets are available for determination of the 
overall PSD but the RXTE monitoring observations provide 
the only information on timescales longer than ~days. We 
therefore begin by fitting a simple power law to the PSD 
from the combined 6-hr, 2-d and total RXTE lightcurves. We 
retain the intrinsic 6-hr and 2-d resolution for the first two 
lightcurves but bin the total lightcurve up to 28-d resolution. 

The resulting PSD is well fitted (fit probability = 68 
percent) by a slope of 0.9±0.15 (90 perc ent confidence inter - 
vals are used throughout this paper). In llJttlev et al.l ( 2002! ) 
we describe how we determine the goodness of fit from the 
simulations. A fit probability of 68 percent means that the 
model is rejected at 32 percent confidence. The fit is shown 
in Fig. |5] For a GBH in the low-hard state, we do not ex¬ 
pect the region of slope ~ — 1 to extend for more than one or 
1.5 decades, before flattening to a slope near zero. Although 
this fit would not be very sensitive to breaks near to the end 
of the frequency spectral range covered, nonetheless the fit 
to a simple power law, of slope close to —1 over approxi¬ 
mately three decades, from ~ 1CP 8 to ~ 10~ 5 Hz, suggests 
that MCG-6-30-15 is the analogue of an XRB in a high-soft 
state rather than in a low-hard state. 


In order to determine properly the overall long and short 
timescale PSD, and hence measure any break frequency and 
high frequency PSD slope, and refine measurements of the 
low frequency PSD slope, we must combine the RXTE mon¬ 
itoring observations with observations which sample shorter 
timescales. 

Both the RXTE long looks and the XMM-Newton long 
observations provide information on shorter timescales. The 
RXTE long looks provide a good determination of the PSD 
on timescales shorter than the length of the individual seg¬ 
ments (~ 3000s). However on timescales between ~ 3000s 
and Mew days, the PSD is distorted by the many gaps in 
the datasets. Our PSRESP analysis software is able to cope 
with those gaps, but they do give the software considerable 
work to do as the dirty PSD differs a good deal from any 
likely underlying true PSD. Thus the errors are larger than 
if there are fewer gaps. The XMM-Newton data are contin¬ 
uous and so the XMM-Newton PSD suffers negligibly from 
distortion. However the sensitivity of XMM-Newton is less 
than that of RXTE in the 2-10 keV range, which is the only 
energy band in which the long timescale PSD is determined. 
Nonetheless, the long XMM-Newton observations of MCG- 
6-30-15, which total approximately three times as long as 
those of NGC4051, do allow a reasonable determination of 
the high frequency PSD. 

As discussed in lM c Hardv et al . (12004 1. the mean photon 
energy of the 4-10 keV XMM-Newton band is approximately 
the same as that of the RXTE 2-10 keV band. In order 
to derive a PSD of uniform energy at all frequencies we 
have therefore included the 4-10 keV XMM-Newton PSD 
with the RXTE monitoring datasets described above in our 
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PSD simulation process. The high-soft state PSD of Cyg X- 
1 is better described by a bending powerlaw model than a 
sharply broken powerlaw model, and the PSD of NGC4051 
is also slightly better fit by the bending powerlaw model 



where q_l and oh are the low and high frequency power- 
law slopes respectively and Ub is the bend, or break, fre¬ 
quency. As stated in the Introduction, here we use the con¬ 
vention that a slope decreasing towards high frequencies will 
be defined by P(u ) oc u a where a is a negative number, eg 
a = - 2 . 1 

We fit this model to the combined RXTE and XMM- 
Newton observations and the result is shown in Figs. © and 
□ In Fig. ©we plot the log of power vs. the log of frequency, 
as in Fig.© In both of these figures the observed dirty PSDs 
from the various constituent lightcurves are given by the 
jagged lines. The model, distorted by the effects of sampling, 
red noise leak and aliasing, is given by the points with error- 
bars, and the underlying, undistorted, model is given by the 
smooth dashed line. In Fig. ©we unfold the observed dirty 
PSD from the distorting effect of the sampling in order to 
produce the closest approximation that we can to the true 
underlying PSD. Thus displacements of the observed PSD 
from the model distorted PSD are translated into displace¬ 
ments from the best-fit underlying model PSD. The tech¬ 
nique is identical to that used to deconvolve energy spectra 
from count rate (ie instrumental) spectra in standard X- 
ray spectral fitting. In this case the error is translated onto 
the observed datapoints and the underlying best-fit model 
is shown as a continuous line. As in standard X-ray spec¬ 
tral fitting, the position of the datapoints does depend on 
the shape of the assumed, best-fitting, model. Note that in 
Fig.© we plot frequency x power so that a horizontal line 
would represent Power (u) oc W 1 and equal power in each 
decade. Similar plots f or NGC4051, NG C3516 and Cyg X-l 
are shown in Fig. 18 of lM c Hardv' et al . (1200411 . 

The combined PSD is well fitted (P=45 percent) by 
this model with low frequency slope a l = —0.8lg'j 6 , high 
frequency slope aH = — 1.98jtg'4Q an d break frequency 
ub = 6.OlJ 0 x 10~ 5 Hz. A sharply broken powerlaw, although 
tolerable, is not such a good fit (P=14 percent) but the best 
fit parameters are similar to those of the smoothly bending 
model (q_l = —0.8, olh = —1.74 and vb = 3.2 x 10~ 5 Hz). 
As is often the case, the break frequency is slightly lower for 
the sharply bending model than for the smoothly bending 
model but, in this case, the errors are not well determined. 
(Note that our software performs a simple grid search so 
the values given here are those of the minimum point on 
the grid and so may differ very slightly from the optimum 
values which may lie slightly off any particular grid point.) 
The confidence contours for these parameters are shown in 
Figs. © © and run 

1 Th is convent i on is used throughout the Tables and text of 
|M C Hardv et alJ i2004T) . No te that when we defin ed the formula 
for a bending powerlaw in lM c Hardv et H . 6553) (top of second 
column, page 788), we forgot to make the formula consistent with 
the rest of the text and Tables and so, only in that one for¬ 
mula, a slope decreasing towards high frequencies is defined by 
P(y) oc v~ a where a. is a positive number, eg a = +2. 



Figure 6. Combined 2-10 keV RXTE and 4-10 keV XMAd- 
Newton PSD of MCG-6-30-15. The lines and datapoints are as 
described in Fig. \S\ 



Figure 7. Unfolded combined 2-10 keV RXTE and 4-10 keV 
XMM-Newton PSD of MCG-6-30-15. Here the errors have been 
translated onto the datapoints and the continous line represents 
the best-fit underlying PSD. Note that here we plot frequency 
xPower. See text for details. 

The break frequency for ,, assumi ng a high state, bend - 
ing powerlaw model, is 22.9 ± 1.5Hz \M C Hardv et ~aIJ l2004ll 
(or 13.9±0.8Hz for a sharply breaking high state model). As¬ 
suming linear scaling of black hole mass wi th fr eque ncy, and 
a black hole mass of lOMgfor a llHerrero et al . ljU), we de¬ 
rive a black hole mass for MCG-6-30-15 of 3.61J, 8 x 10 6 Mq. 

We note that, although con sistent w ithin th e error s with 
the break frequency derived by I Vaughan et all (1200 30 using 
purely XMM-Newton observations, our break frequency is 
slightly lower. The reason for the difference is that Vaughan 
et al. were unable to measure the PSD slope below the break 
and so assumed a value of -1. However our data show that 
the slope is actually slightly flatter, which leads to a lower 
break frequency. 
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Figure 8. 68 percent, 90 percent and 99 percent confidence con¬ 
tours for high frequency slope, olh-, and break frequency, i/jg, for 
bending powerlaw fit to the combined RXTE and XMM-Newton 
4-10 keV PSD shown in Fig. EJ Note we plot —olh- 


Figure 10. 68 percent, 90 percent and 99 percent confidence 
contours for low and high frequency slopes, oll and olh respec¬ 
tively for bending powerlaw fit to the combined RXTE and XMM- 
Newton 4-10 keV PSD shown in Fig.|b] Note we again plot —oll 
and —olh- 



Figure 9. 68 percent, 90 percent and 99 percent confidence con¬ 
tours for low frequency slope, oll , and break frequency, vb, for 
bending powerlaw fit to the combined RXTE and XMM-Newton 
4-10 keV PSD shown in Fig.|b] Note we plot —oll 

3.2.2 RXTE monitoring observations and low energy 
XMM-Newton observations 

In lM c Hardv et alJ d2004ll we refined our determination of 
the break frequency in NGC4051 by combining RXTE mon¬ 
itoring observations with continuous XMM-Newton obser¬ 
vations in the 0.1-2 keV band, where the break is better 
defined. We assumed that the break frequency was indepen¬ 
dent of energy. We have carried out the same procedure here 
(again taking proper account of the different PSD normali¬ 
sations between the 0.1-2 and 2-10 keV bands) and find that 
a L = -0.8+°;,, OLH = — 2 . 5^4 and v B = 7 . 6 + 3 ° x 10“ 5 Hz. 
The fit probability is 67 percent. Apart from the steeper 
high energy slope which is well known at low energies 


llVaughan et alJ |2003t lM c Hardv et al.ll2004lf . the other fit 
parameters are very similar to those obtained when using 
the RXTE and XMM-Newton 4-10 keV observations. Thus 
within the errors, we find no evidence for a change of break 
frequency with photon energy. The implied black hole mass 
is more tightly constrained by the smaller errors on the break 
frequency and is 2.9+,'g x 10 6 A7 q. 

We fitted a sharply breaking powerlaw model to the 
combined RXTE and XMM-Newton 0.1-2 keV data. As with 
the combined RXTE and 4-10 keV XMM-Newton data, the 
fit parameters are almost exactly the same as for the bend¬ 
ing powerlaw model, but the fit probability is worse (19 per¬ 
cent). The fact that a bending powerlaw is a better fit than 
a sharply bending powerlaw to the PSD of MCG-6-30-15, 
as it is to the PSD of in its high state, strengthens our con¬ 
clusion that MCG-6-30-15 is the analogue of a GBH in the 
high, rather than low, state. 

3.2.3 RXTE monitoring observations and RXTE long 
looks 

We have also determined the overall PSD shape by includ¬ 
ing the RXTE long looks binned up to 128s time resolution 
with the RXTE monitoring observations, and also includ¬ 
ing a very high frequency PSD (> 10~ 3 Hz) made from the 
individual segments of the RXTE long looks. The fit param¬ 
eters are approximately the same as in the combined RXTE 
and XMM-Newton fit ( oil = — O. 8+0151 a H = —2.1+g g and 
v B = 5 . 0^4 x 10~ 5 Hz). The fit is formally good (P=74 per¬ 
cent) but the errors are rather large because of the diffi¬ 
culty of coping with the many gaps in the RXTE long look 
lightcurves. 

3.3 Search for a second, lower frequency, break 

We searched for a second, lower frequency break, using the 
RXTE and XMM-Newton 4-10 keV data. We assumed a low 
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L°9 ^high break 

Figure 11. Confidence contours for the high and possible low 
break frequencies. The 68 percent contour is solid, the 90 per¬ 
cent contour is dashed and the 99 percent contour is dot-dashed. 
The straight dashed lines, with labels of 10 and 100 respectively, 
indicate those ratios between the high and possible low break fre¬ 
quencies. The fits assume a break below the low frequency break 
of 0 and an intermediate slope between the low and high frequency 
breaks of 0.8. Both break frequencies were allowed to vary, as was 
the slope above the upper break, which was measured at 1.95T0.3. 

state PSD model, fixing the slope below the lower break at 
0 , and the slope between the lower and upper breaks at - 
0 .8, i.e. the value we measure in the same frequency range 
assuming a high state model. We allowed the two break fre¬ 
quencies and the slope above the upper break to vary. The 
best fit frequency for the higher break and for the slope 
above the higher break was the same as in our high state fit 
and the best fit for the lower break frequency (6 x 10 -9 Hz) 
was almost at the end of the fitting range. The fit probabil¬ 
ity was 43 percent. In Fig. mi we show the 68 percent, 90 
percent and 99 percent confidence ranges for the two break 
frequencies. As the typical ratio of the two break frequencies 
for a low state GBH is between 10 and 30, we plot, as dashed 
straight lines, the limits corresponding to a ratios of 10 and 
of 100. The best fit ratio is over 1000, but a ratio of 100 is 
just within the edge of the 90 percent confidence region. 

We repeated the above analysis after fixing the interme¬ 
diate slope at —1.0 and obtained almost identical frequency 
ratios although the high break frequency was then about 
half a decade higher in frequency. 

On the basis of its PSD, the above analysis strongly 
suggest that MCG-6-30-15 is the analogue of a high state 
GBH, although a low state cannot be entirely ruled out. 


4 MASS DETERMINATION: ABSORPTION 
LINE VELOCITY DISPERSION 

As discussed in the Introduction, the mass of the black hole 
in MCG-6-30-15 is not well determined. In particular there 
have been no reverberation mapping observations, and no 
measurement of the stellar velocity dispersion, the two tech¬ 
niques widely regarded as giving the most reliable measure¬ 


ment of black hole masses. However a reliable measurement 
of the mass in MCG-6-30-15 is important for any discussion 
of mass/timescale scalings in AGN and for any discussions 
of whether AGN are in low or high states. In this section we 
therefore present stellar velocity dispersion measurements 
from which we estimate the black hole mass. In subsequent 
sections (Sec 0 and Sec^Jl we present secondary determina¬ 
tions of the black hole mass, based on the width of the H (3 
emission line and on photoionisation calculations. We find 
that all these optically based mass determination methods 
give consistent answers. 

4.1 Observations 

The observations of MCG-6-30-15 described here were taken 
by the service programme of the 3.6m Anglo-Australian 
Telescope on 2002 June 5. The RGO Spectrograph was 
used, with the 25cm camera plus the 1200R grating (blazed 
at 7500 A). The detector is an EEV2 CCD, windowed to 
600 x 4096 pixels for the science observations and 150 x 4096 
pixels for the standard stars. A slit width of 1" (0.15mm) 
was chosen, resulting in a spectrum of resolution R ~ 9500 
covering a wavelength range of 1000 A. The slit was placed 
along the major axis of the galaxy, at a position angle of 
116°. 

The measured widths of the Call triplet lines are used 
to estimate the black hole mass. At the redshift of MCG-6- 
30-15, z = 0.00775, the Can lines are shifted to AA8566, 8610 
and 8731 A, and therefore the central wavelength was cho¬ 
sen to be 8620 A. The total wavelength coverage of 1000 A, 
spanning the wavelength range 8120 — 9120 A, allows accu¬ 
rate determination of the continuum on either side of the 
lines, and therefore correct measurement of the line widths. 

Observations of late-type giant stars, which have very 
narrow lines, and which dominate the stellar light of nearby 
galaxy bulges, were performed in order to measure the in¬ 
strumental broadening and provide templates for the cross¬ 
correlation analysis. The K0 ill standard stars HD 117927 
and HD 118131 were observed, with the same set-up and 
wavelength coverage as MCG-6-30-15, since the rest wave¬ 
length of the Call triplet (AA8498, 8542, 8662A) lies within 
the wavelength range covered. 

Dedicated flat field exposures were taken before and 
after the science and standard star observations, for accurate 
de-fringing, and argon and neon arc spectra were taken for 
the wavelength calibration. 

4.2 Data Reduction 

The raw data were reduced with IRAF using standard bias 
subtraction, flat fielding and cosmic ray removal techniques. 
At the red wavelengths observed here, the EEV2 CCD suf¬ 
fers from fringing at a level of ~ 5 percent. A comparison of 
the extracted science spectra with equivalent spectra taken 
from the flat field frames showed that the fringes were effi¬ 
ciently removed from the science frames, and did not con¬ 
tribute any spurious features. Curvature in the spatial direc¬ 
tion was rectified using bright night sky lines in the object 
frames. The narrow window used for the standard star ob¬ 
servations meant that in this case, rectification was unneces¬ 
sary. The continua of MCG-6-30-15 and the standard stars 
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Figure 12. Normalised spectra of MCG-6-30-15 (top panel) and 
a template star (bottom), with the position of the calcium triplet 
feature and broad Ol A8446 emission line marked. 

were sufficiently bright that curvature of the spectra in the 
dispersion direction could be traced and rectified. The total 
exposure time for MCG-6-30-15 was 1 hr, broken into two 
30min integrations. Wavelength calibration was performed 
using the arc exposures bracketing each science frame. This 
calibration was checked by measuring the resulting wave¬ 
length of night sky emission lines, and no further correction 
was required. 

The nuclear spectrum of MCG-6-30-15 was extracted 
from the central 5.6 pixels of the trace, which corresponds 
to a 2.4" x 1" aper ture. To be compatible with the Mbh — cr* 
relation defined bv lFerrarese fc MerritfJ feoQ Q). the length of 
this aperture was chosen to correspond to r e /8, where r e is 
the effective radius of the galaxy, measured from the profile 
of the spectral data, compressed along the dispersion direc¬ 
tion. We note however that the choice of extraction aperture 
has only a very small effect on the results obtained (Ho, pri¬ 
vate communication). Sky subtraction was performed using 
regions either side of the galaxy, 50 pixels (21.5") wide, at a 
distance of 50 pixels (21.5") from the centre of the galaxy. 
The resulting spectrum is shown in Fig. EH normalised by 
a spline fit to the continuum to emphasize the emission and 
absorption features. 

In a similar analysis, iFilippenko fc lid J 200. 'ft) detect 
Paschen emission lines in the spectrum of NGC4395, adja¬ 
cent to each of the calcium absorption features. For MCG- 
6-30-15, there is perhaps a suggestion of a broad Paschen 15 
line at ~ 8615A, but it is not conclusive. Even if present, 
the strength of such a line is very small compared to those 
detected in NGC4395, and as such, would not significantly 
affect our results. 


4.3 Data Analysis 

The calcium absorption features in a galaxy spectrum are 
broadened compared to those of an individual star due to 
the bulk motion of the stars in the galaxy. The extent of 
this Doppler broadening is related to the mass of the central 
black hole, giving rise to the observed tight relationship be¬ 
tween the two quantities. The region of the galaxy spectrum 
of MCG-6-30-15 containing the calcium triplet is therefore 


cross-correlated with the equivalent region of the KOlll stel¬ 
lar templates using the IRAF task fxcor. The task fits a 
spline function to the continuum of each input spectrum, 
and cross-correlates the resulting features, giving as output 
a velocity measurement, due to the redshift of MCG-6-30- 
15, and a velocity width, due to the Doppler effect of the 
motion of stars in the galaxy caused by the gravitational 
influence of the black hole. 

The width of the cross-correlation peak function cal¬ 
culated by fxcor measures the combined Doppler and in¬ 
strumental broadening, therefore in order to ascertain the 
Doppler broadening alone, the stellar templates themselves 
were broadened using a Gaussian of various widths, and 
cross-correlated with the unbroadened templates. The width 
of the Gaussian which gave rise to a cross-correlation peak 
width equal to that of MCG-6-30-15 provides our mea¬ 
surement of the true velocity dispersion of the galaxy. We 
found that broadening the stellar templates by 11 pixels 
gave the best match to the cross-correlation peak width 
of MCG-6-30-15. The velocity dispersion of the spectra is 
8.5kms -1 pix -1 , giving a velocity dispersion for MCG-6- 
30-15 of a = 93.5kins _1 . 

The statistical errors on this measurement are very 
small, and therefore d o not provide a realist ic estimate of 
the true errors involved. iFerrarese et ali (; 200ll) estimate that 
the systematic uncertainties involved in making such mea¬ 
surements are of order 15 percent. We therefore undertook 
a number of tests in order to ascertain the likely error in our 
result. First, since we have three standard star observations 
(two of HD 117927 and one of HD 118131), each spectrum 
was broadened in turn and cross-correlated with the two un¬ 
broadened stellar templates to estimate the variation due to 
the use of differing standards. The results usually varied by 
less than one percent, and never more than two percent. Sec¬ 
ondly, to assess the possible impact of the proximity of the 
Ol line to the first feature in the calcium triplet, we restrict 
the region of the spectra used in the cross-correlation analy¬ 
sis to that containing only two of the three calcium absorp¬ 
tion lines. By changing the wavelength range over which the 
cross-correlation analysis is performed, a maximum varia¬ 
tion of ±4 percent in the peak width was measured. Thirdly, 
we vary the amount of smoothing of the stellar templates 
required to match the result from MCG-6-30-15, and find 
that all the results of the above tests are easily bracketed 
by taking a Gaussian width of 11 ± 1 pixels, which is equiv¬ 
alent to cr = 93.5 ± 8.5 kms^ 1 (±9.1 percent). We therefore 
take this as our conservative estimate of the error in our 
measurement, and note that this is small compared to the 
uncertainty in the relationship between this quantity and 
the inferred black hole mass, as described in the following 
section. 

4.4 Black Hole Mass Estimate 

The largest source of error in estimating the mass of the 
black hole from the width of the absorption lines is not 
uncertainty in the measurement of the width but uncer¬ 
tainty in the Mbh — cr relation itself. The two main groups 
concerned with deriving the relationship use slightly dif¬ 
ferent observational methods and fit slightly d ifferen t re¬ 
lationships to the resultant data (see iMerritt fo Ferra.resd 
fiooit iTremaine et~aHf200 2i for full discussions). Recently 
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iGreene et al. ( 12004 ) have shown that AGN with very low 
black hole masses (down to ~ 1O 4 M 0 ) fit reasona bly well 
onto the relationship given hv lTremaine et al ] (120021) . i.e., 


«»» - 1.35iS:S x 10 6 Mq 


4.02(±0.32) 


Our observations then impl y Mbh = 6.3 l 3 ° x 1O 6 M 0 . As 
can be seen from Fig. 4 of lGreene et al. ( (2004 1. the spread 
in the datapoints implies an uncertainty in the derived mass 
of at least 50 per cent. 

An alternative version o f this relationship has been 
given bv lMerritt fc Ferrarese 1200 11. Using the most r ecent 
version of this alternative relationship llFerrareseil2002 l i.e. 


/ a x4.58(±0.52) 

Mbh = 1.66(±0.32) x 1O 8 M 0 —- 

V 200 km s / 

we derive Mbh = S-ljl®'® x 1O 6 M 0 , which is entirely con- 
sistent wit h th e mas s derived using the relationship of 
iTremaine et alJ J200 2). 



Figure 13. HST/STIS spectrum of MCG-6-30-15, showing 
broad and narrow H/3 components, and [Om] AA 4959,5007A 
emission lines. 


5 MASS DETERMINATION: EMISSION LINE 
WIDTH 


An alternative estimate for Mbh can be o btai ned us ing the 
empirical relationships found by iKasoi et all ( 2000) in re¬ 
verberation studies of nearby Seyfert 1 galaxies and quasars 
between Mbh, the velocity dispersion of the broad emission¬ 
line gas, Vfwhm, 2 and the effective size, Rblr, of the broad 
line emitting region, i.e.: 


Mbh = 1.464 x 10* 


( Rblr \ 

( Vfwhm ) 

V It — days J 

V 10 3 km s _1 / 


M 0 . 


Rblr is determined from the measured delay between the 
continuum and emission-line variations and is related to the 
galaxy subtracted continuum luminosity at A5100A by 


Rblr = 32.9)) 2 9 


/ AL a (5100A) 
Y 10 44 erg s _1 


0.7 


[Reynolds et all (ll997ll discuss in detail the optical spectrum 
of MCG-6-30-15 and, in their Fig. 4, derive an estimate 
of the intrinsic non-stellar optical flux (ie uF u , or A F\) at 
5100Aof ~6x 10 -11 ergs cm 2 s 1 . Assuming the standard 
flat cosmology of Ho = 75 km s” 1 Mpc -1 and A matter = 
0.3, with 2 = 0.007749, this flux equates to a rest-frame 
luminosity of \L\ = 7.2 x 10 42 erg s _1 and yields Rblr = 
5.2 lt-day. 

We note that there is considerable scatter in_the re¬ 
lation ship between Rblr and A F\ at MOOA^JVes^reaardl 
(2002) repeated the empirical study of|Kasg^MilJj|j^^7 
using the linear regression techniques of tAkritas fc Bershadvl 
( 1996 ) and obtained a slightly different relationship. How¬ 
ever the derived value of Rblr, ie 5.36 lt-day is, within the 
er rors, identica l to th e value derived usin g the relationship 
of lKasoi et alJ ll200fltL IVestergaardl ll2002l) derives the rela¬ 
tionship using only broad line Seyfert 1 galaxies. The narrow 


2 While one might expect Vfwhm of the root mean square pro¬ 
file to better represent the velocity dispersion of the variable 
part of the BLR, and thus show a better correspondence with 
Rbir. in practice there is little difference in the derived masses 
iKaspi et al lToOOlb 


line Seyfert galaxy, NGC4051, is an outlier and i s not used in 
d erivi ng the relationship, although it is used bv lKasoi et alJ 
( 2000 ). 

Using archival HST/STIS data for MCG-6-30-15 
(Fig. uni we have measured (rest frame) Vfwhm = 32.9 ± 
1.9A (= 2020 ± 120 km s _1 ) which, together with Rblr, 
gives a virial mass for the black hole in MCG-6-30-15 of 
3 x 1O 6 M 0 . 


6 MASS DETERMINATION: 

PHOTOIONISATION 

The size of the BLR may also be found from photoioniza¬ 
tion calculations, but with relatively large uncertainty. The 
calculated BLR size, combined with the velocity at FWHM 
of the H (3 line profile can then be used to provide an esti¬ 
mate for the virial mass of the central black hole. From a 
sam ple of 17 Seyfert 1 galaxies and 2 quasars. IWandel et alJ 
( 19991) found an approximately linear relationship between 
photoionization mass estimates and reverberation mass es¬ 
timates, suggesting that photoionization calculations might 
providing a route for mass determinations in systems with 
only single epoch s pectra l observations. The empirical rela¬ 
tionship derived bv lWandel et alJ f 19991) is: 


Mbh = 2.8 x 10 6 / 


/ L44 \ / Vfwhm \ 2 

\Unio) VKPkins” 1 / 


M 0 . 


Here U is the ionization parameter (the dimensionless ratio 
of photon to gas density), nio is the electron density in units 

of 10 10 cm -3 , and / is the product fkfj/ 2 E ^ , where fk 
is a factor relating the effective velocity dispersion to the 
projected velocity dispersion, relates the observed lumi¬ 
nosity L to the ionizing luminosity, and E = E/1 ryd. If we 
assume a BLR ionization parameter and gas density typical 
for nearby Seyfert 1 gal axies (U = 0.1, nio = 10), and a 
weighted /-value ~ 1.45 IWandel et all (1999) we derive a 
black hole mass of 4.5 x 1O 6 M 0 . 
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7 DISCUSSION 

7.1 Summary of Mass Determinations For 
MCG-6-30-15 

Although considerable systematic uncertainties and large as¬ 
sumptions (eg in nio an d U) are involved in their derivations, 
the various optical measurements of the black hole mass in 
MCG-6-30-15, including the revised estimate based on the 
bulge mass, are in reasonable agreement. All values he be¬ 
tween ~ 3 and ~ 6 x 10 6 Mq. We cannot tell which mea¬ 
surement is correct and so, as a working value, we take the 
middle of the range, ie ~ 4.5 x 10 6 MQand adopt an error 
equal to the spread in the measurements (ie 3 x 10 6 Mq). We 
note that although there are again large uncertainties in the 
mass derived from the PSD, the most tightly constrained 
mass, ie that derived from a combination of RXTE and low 
energy XMM-Newton observations and assuming linear scal¬ 
ing of break timescale with mass from Cyg X-l in the high 
state, ie ~ 2.9^'® x iV Mq, is in good agreement with the 
mass as determined by optical methods. 

The average 2-10 keV X-ray flux of MCG-6-30-15 from 
our RXTE monitoring is 5.9 x 10~ n ergs cm -2 s~ 4 , 
which corresponds to a luminosity of ~ 7 x 10 42 ergs 
s~ 4 . For an assume d X-ra y/bo l ometric co rrection of 27 
llPadovani fc Rafanellilll98Sl : Elvis et al.ll994ll . the bolomet¬ 
ric luminosity is 1.9 x 10 44 ergs s _1 , implying that, for a mass 
of ~ 4.5x 10 6 Mq, MCG-6-30-15 is radiating at ~ 40 percent 
of its Eddington luminosity. For narrow line Seyfert galax¬ 
ies, the bolometric correction should probably be less than 
27 but probably still greater than 10, so MCG-6-30-15 is 
almost certainly radiating at > 10 percent of its Eddington 
luminosity. A low state interpretation of the PSD which, al¬ 
though very unlikely, cannot be ruled out entirely, implies 
a black hole mass of ~ 5 x 1O 5 M0, thus requiring a super- 
Eddington luminosity. 


7.2 X-ray Selected AGN and High State PSDs 

Using our newly derived black hole mass, and slightly refined 
break timescale, we plot (Fig. II41 MCG-6-30-15 on a re¬ 
vised version of the break timescale/black hole mass diagram 
which we presented in M c Hardy et al. (2004). Although 
bending, rather than sharply breaking, powerlaws fit the 
PSDs best, we still use timescales derived from sharp breaks 
in this diagram as we do not yet have timescales derived from 
bending powerlaw fits to all of the AGN. Like NGC4051, 
MCG-6-30-15 also sits above the high state lin e. T h e break 
timesc ales are taken from the compilation in|MMIcUx^^^iI 
|222^) with the addition of NGC4395_from^^^^E^^^L 
12QQa), with mass es timate from iFilippenko fc Hoi <20031. 
iKraemer et alJ lll999l) present UV and optical spectral of 
NGC4395 and show broad, although rather weak, wings to 
H/3. They estimate FWHM^1500 km s _1 , but do not give 
an error. We are therefore unsure whether to class it as a 
broad or narrow line Seyfert galaxy and so mark it by an 
open triangle in Fig. 1141 We also include NGC3227 from Utt- 
ley and M c Hardy (in preparation). Uttley and M c Hardy also 
discuss NGC5506 in considerable detail but we do not in¬ 
clude it here as its mass is highly uncertain. 

Where available (ie Fairall9, NGC3227, NGC3783, 
NGC4051, NGC4151, NGC5548) we use reverberation 


masses from the compilation of lPeterson et al~ (2004 .1. No 
reverberation masses are available for the other AGN, ie the 
narrow line Seyfert Is, and so masses are derived from stel¬ 
lar velocity dispersion measurements. For some objects al¬ 
ternative masses are available from the emission line width 
(eg 6.3 x 10 s Mq for Mkn766 from Botte et al. 2005) but 
to reduce the number of variables we restrict ourselves to 
the velocity dispersion mass (3.5 x 10 6 Mq for Mkn766, us¬ 
ing the calibration of Tremaine et al. 2002). C urrent work 
(e.g. lFerrares~ et alJl200 L; IPeterson et alj|2004h shows that 
masses derived from reverberation mapping and velocity dis¬ 
persion are consistent. An exception is Akn564. As in some 
other NLSls, t h e Ca n trip let fines in Akn564 are in emis¬ 
sion ( Ivan Groningenlll99 3). rather than absorption, and so 
cannot be used to determine the black hole mass. In this 
case the width of the [O III] 5007 emission line is used as 
a substitute for the width of the stellar absorptio n l ines in 
order to estimate the black hole mass dBotte et alJ 12004 1. 
The width of the [O III] 5007 emi ssion line does correlate, 
although with considerable scatter dBorosonll2003l) , with the 
width of the Cal l trip l e ste llar absorption lines. However it 
has been noted dBotte et al.l200^i that the [O III] lines tend 
to be wider than the Ca absorption lines. Therefore, assum¬ 
ing that the Ca absorption lines represent the black hole 
mass more accurately, the [O III] fine will typically give an 
overestimate of the black hole mass, hence the upper limit 
on Akn564 in Fig. ITT! 

Due to recalibration to better fit the Mbh — o* re- 
lationship, masses in the sample of IPeterson et aid ('2 004fl 
ha ve generally inc reased from the earlier estimates used 
in Im Ha.rdv et alJ ( 20041 ). Also, the b lack hole mass esti¬ 
mate for NGC4051 fromJPeterson_et ah|j |2004h is higher 
th an the estimate f romls herTune^^m2003il which we used 
in lM c Hardv et all (l2004h . Here, for consist e ncy, w e use the 
black hole mass values from IPetersoi]_j4 i ja : lj _(|2004)jwherever 
available. Thus all AGN from tho lPeulmon et all (l2004h sam¬ 
ple have moved further away from the low-state line and to¬ 
wards the high state line. With the exception of NGC4151, 
all AGN now lie above the low state line. Although the 
‘broad line’ Seyfert galaxies are, in general, closer to the 
low state line than the ‘narrow line’ Seyfert Is, it is possible 
that all AGN considered here might be high state systems. 
Thus although rh may well have an important effect on de¬ 
termining the break timescale, the accretion rate in all AGN 
studied here ma y be high enough to m ake them ‘high state’ 
systems. Indeed IPeterson et all •: 2(1(1 ill show that the aver¬ 
age accretion rate for the current sample of AGN is roughly 
10 percent of the Eddington rate, which exceeds the ~ 2 per¬ 
cent rate which typifies the tr ansiti o n to the high state in 
galactic X-ray binary systems dMaccarond 12003 1. For com¬ 
pleteness we also include the break timescale of NGC3227 
from Uttley and M c Hardy (2005, in preparation). NGC3227 
is an interesting galaxy, having broad permitted lines, a hard 
X-ray spectrum, and a high state PSD. It is discussed in de¬ 
tail by Uttley and M c Hardy (2005, in preparation). 

An additional indicator of the ‘state’ of an accreting 
black ho le is given by the ratio of its X-ray to radio luminos¬ 
ity fe.g. iGa.llo et al.ll2003l : lFenderll2fifnl 1200,‘ji . For Galactic 
black hole systems, a high X-ray/radio ratio signifies a high 
accretion rate and a ‘high’ state system. A low X-ray/radio 
ratio signifies a low accretion rate and a jet-dominated ‘low’ 
state system. The core radio flux of MCG-6-30-15 is very 
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low, ~ lmjy at 5GHz (llJlvestad & Wilsonlll984ll . and very 
similar to that of NGC4051 (M c Hardyet al, in preparation). 
The X-ray fluxes, and hence the X-ray/radio ratios, of both 
galaxies are high, again indicating ‘high’ state systems. 

7.2.1 A Possible Selection Effect 

It is interesting to note that the few AGN with sufficiently 
good PSDs to distinguish between high and low states 
(NGC4051, MCG-6-30-15, NGC3227), all have high state 
PSDs . The AGN which we, and others le.g. lMarkowitz et a.lJ 
2003), have been monitoring are mainly taken from the X- 
ray bright AGN which are vi sible in the bright-flu x limit, 
all-sky, X-ray catalogues (e.s. iMcHardv et all 1981 1. As X- 
ray flux rather than, eg, radio flux, is the only selection 
criterion in such surveys, there is a strong selection effect 
towards selecting high state AGN. 

AGN with low state X-ray PSDs are likely to be found 
amongst those A GN which have rela tively more luminous ra¬ 
dio emission (cf lMerloni et al.ll2003tl . As X-ray emission will 
not be the only parameter on which we select such AGN, 
we may expect them to be, in general, fainter in X-rays 
than those studied presently. Such AGN should be present 
in, eg, the sample of radio and X-ray bright objects se¬ 
lected from the ROSAT and VLA FIRST all-sky catalogues 
iferinkmann et aDl200f)l) . 

The ob servations from RXTE presented here, and else - 
where (e.g. iMarkowitz et, alJ l2003t lM c Hardv et alJ [200411 . 
demonstrate the tremendous importance of long timescale 
(years/decades) X-ray monitoring observations for our un¬ 
derstanding of AGN. RXTE has revolutionised our under¬ 
standing and opened up many new exciting areas of study, 
and hopefully it will continue to operate for many years to 
come. However, in the longer term, it is crucial that a sen¬ 
sitive (^mCrab in a few hours) all sky X-ray monitor, with 
a long lifetime (^decade) is launched to carry on this excit¬ 
ing work, otherwise this newly emerging field will die with 
RXTE . 
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Figure 14. Black hole mass vs. PSD Break Timescale. Broad line 
Seyfert galaxies are shown as filled circles, narrow line Seyfert 
galaxies are shown as open circles and NGC4395 is shown as 
an open triangle. The dashed line to the right represents linear 
scaling of break timescale with mass from Cyg X-l in its low state. 
The dot-dash line, on the left, represents linear scaling of break 
timescale with mass from Cyg X-l in its high state. NGC4395 is 
plotted as an open triangle (see text). The arrow labelled with m 
indicates the way in which the relationship of break timescale vs. 
mass might move with increasing accretion rate. 
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